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Abstract. Consider a particle moving through a random medium, which consists of spher- 
ical obstacles, randomly distributed in R d . The particle is accelerated by a constant external 
field; when colliding with an obstacle, the particle inelastically reflects. We study the asymp- 
totics of X(t), which denotes the position of the particle at time t, as t — > oo. The result is 
a functional limit theorem for X(t). 
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1. The Lorentz model and the problem 

1.1. The motivation. Consider a spherical particle moving in a random medium. The 
medium consists of immobile spherical obstacles of equal radii, randomly distributed in IR 3 . 
The particle is accelerated by an external field providing constant acceleration a. At a 
collision with an obstacle, the particle's speed v changes to 

v — (1 + a)(v , v)v, 

where < a < 1 is the restitution coefficient and v is the inner unit normal to the obstacle 
at the point of collision. 




an obstacle 
the particle 



This mapping changes only the normal component of v, i.e., (v, u)u, which is multiplied by 
—a. On the figure, the dotted line indicates the trajectory of the particle's center. 

This model is often named after Hendrik Lorentz who introduced it (see [8]) in order to 
describe conductivity in metals. Lorentz studied the case of elastic collisions, with a = 1; 
the generalization for non-elastic collisions could be found, e.g., in |15j . 

In physics the Lorentz model is used to describe the motion of a particle in a medium if 
the particle's mass is negligible with respect to masses of the medium's particles (obstacles). 
Indeed, in this case we can assume that the obstacles have infinite masses and thus they 
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remain immobile at collisions. For example, the model in question describes well the motion 
of electrons in helium, see [16J. 

In case of elastic collisions, the Lorentz model is a billiard type model; recall that a 
billiard is the following model of motion of a particle in an arbitrary region with smooth 
boundary: at a collision with the boundary, the particle elastically reflects. Indeed, we can 
consider the region whose boundary consists of the obstacles' boundaries. The case when 
the obstacles are located periodically, is especially interesting (see [I], [12]); such billiards 
are called Lorentz periodic gases. The main tool for studying billiards is ergodic theory; the 
basic results of this theory and their applications for the Lorentz model could be found in 
[5]. For a detailed review of problems and methods of the billiard theory, see [5] and [T3] . 

Another interesting interpretation of the Lorentz model could be found in [T5] (also, see 
the references therein): a particle percolates through an immobile medium under the constant 
gravity. Considering a big number of such percolating particles and neglecting interactions 
between them, we get a model of mixing of two dry substances, for example, powders. 

There are many physical papers on the Lorentz model. Usually their main goal is to 
study the Boltzmann equation (that is the equation for density p(v,t) of probability that at 
time t the particle's speed is v; this equation could be derived from the law of conservation 
of matter); see pQ, [3], and [9]. For instance, the purpose of [9] is to prove the existence of 
the stationary, i.e., independent of t, solution of the Boltzmann equation in case a < 1. It 
is typical that in [9] the convergence of p(v,t) to the stationary solution, as t — > oo, is not 
discussed. 

1.2. The model. Since the Lorentz model is extremely complicated to analyze, we replace 
it with a simpler one. In doing so, we follow multiple papers on the Lorentz model, for 
example, PQ, [3], [9], and [15]. For the case of elastic collisions, our simplified model coincide 
with one introduced in [11]; in this paper the authors only consider a = 1 and formulate 
their model in a rather different way. We also note that the present simplified model was 
implicitly used in [9]. However, in [9] and [TT] the authors do not discuss how the simpler 
model is derived from the original one. In Section [2] we give such explanation. 

Let us formulate the simplified model. Denote by V n the speed of the particle just 
before the nth collision and denote by r n the random time between the nth and the (n + l)th 
collisions. Let {cr„} n >i C S 2 C M 3 be uniformly distributed unit vectors; let {r] n } n >o be 
exponential random variables with mean A; and let {cr n } n >i, {r] n } n >o be independent. Here 
A > is a parameter signifying the mean free path of the particle. We state that 

1 + oc , N 
V n +i = V n — \y n + \ V n \a n ) + ar n ; (1) 

r n = F(y n - (V n + , Vn) , (2) 

where F : M 3 x R + — ► M is a deterministic function defined as the solution of the equation 

F(v,t) 

\v + as\ds = t. 

In addition, we suppose that at time zero the particle's speed is nonrandom and equals some 
Vq G M 3 . Therefore V\ = v + ar , where r = F(vo,r] Q ) is the random moment of the first 
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collision. We stress that V n form a Markov chain, and the motion of the particle is completely 
defined by this chain. 

Let us briefly compare the new model and the original one. In the new model, at the 
nth collision the speed V n changes to V n — ^Y L (V n + \ V n \a n ). This corresponds to the collision 
with an obstacle, whose inner normal v n directed along the bisectrix of the angle between a n 
and V n (in fact, (V n , z/ n )z/ n = \{V n + |Ki|cr n )). Further, since for every v the function F(v, •) 

is monotone, we have F (v, J* \v + as\ds^ = t. Substituting v := V n — + |V^|cr n ) and 

t := T n , and comparing the resulting equality with (J2]), we get 



Vn 







1 + a 

V n — (V n + \V n \a n ) + as 



ds. 



Thus rj n is the length of the path passed by the particle between the nth and the (n + l)th 
collisions. As the mean free path is A, it is quite natural that rj n are exponential r. vs. with 
mean A. If we know the path's length r\ n and the initial speed V n — ^-^{V n + \ V n \a n ), we can 
find the time r n ; this argument explains ([2]). 

However, there are quite significant disparities between the models. In the new model, 
if the particle collides with an obstacle at some point of space, then the particle does not 
necessarily collide again while coming back to the same point. In other words, after a collision 
the obstacle instantly " disappears" . This happens because of the postulated independence 
of all <J n ,r] n . Thus the medium "changes" in a rather specific way. See more about the 
disparities in Section [2J 

1.3. The problem and the results. We study the asymptotics of X(t), which denotes 
the position of the particle at time t, as t — *■ oo. The case a = 1 was investigated in [11], 
where the authors proved that after a proper normalization the trajectories of X(t) weakly 
converge to a certain diffusion process. 

Let us mention about a very similar model of motion in M 1 : at a collision, the particle's 
speed v always changes to — av. If there is no external field and a = 1, then the times 
between collisions are i.i.d. exponential r.vs., thus X(t) is a well known telegraph process. 
This simplest model of motion is studied in detail in [7] . 

In the current paper, we consider the motion in R 3 , with an external field, and assume 
that a G (0, 1). The purpose of this paper is to prove a functional limit theorem for X(t), 
and thus to sharpen the results of the previous work [2] . We also note that the model in 
question could be easily generalized to obtain the model of motion in M. d . 

Without loss of generality we assume that X(0) = 0. Consider an orthonormal basis of 
M. 3 such that for the acceleration a it is true that a = (0, 0, |a|) T . 
Our main result is the following 

Theorem 1. Suppose < a < 1 and \a\ ^ 0; then there exist constants c\ > and C2, C3 > 
such that for any initial speed v e M 3 , m the space C([0, 1],M 3 ) 

, X(st) — C\ast d I ° 2 !?M x 1 

Yt(s) := V J 1 >¥(■):= c 2 W 2 (-) | , t -> 00, 

V* V c 3 W 3 (-) 
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where W%, W2, and W3 are independent Wiener processes. 

Remark. The constants Ci,c 2; and c 3 depend on the model's parameters a, a, and A. The 
author failed to find these constants in an explicit form. 

Remark. We can easily extend the model defined by the relations flTJ) and (jSJ) to get the 
model of motion in random medium in M. d . Indeed, let a n be uniformly distributed on S^ 1 C 
R d , let a e R d , and let F : R d x R + — > R be defined as above. For this model in R d , 
Theorem^ also holds; the proof is almost the same as for R 3 . For the limit process, we have 
Y = (c 2 Wi, . . . , CzWd-i, c 3 Wd) T , where Wi are independent Wiener processes. 

Moreover, for the motion in R 1 , Theorem^ is valid for the similar model, where at a 
collision, the particle's speed v always changes to —av. 

A nonrigorous explanation of our result could be found in |15j . 

Proving this theorem, we reduce the problem to some statements about a certain Markov 
chain. The most difficult one is that for this chain the functional central limit theorem 
(FCLT) holds. The difficulties arise because the chain has uncountable and noncompact state 
space; moreover, the chain does not satisfy the Doeblin condition, thus the classical results 
(for instance, from [6]) are not applicable. We solve the problem using quite recent results 
(see [10J and references therein), based on stochastic analogues of Lyapunov's functions. 

Therefore, our methods differ significantly from those of [TTJ. Nevertheless, some simi- 
larities could be found on a deep level. Indeed, the proof from [IT] is based on martingale 
theory, while the presented one follows from the FCLT for Markov chains. Recall that the 
FCLT is proved by martingale theory arguments. 

2. The deduction of the simplified Lorentz model 

In this section we consider the original Lorentz model, described in Subsection ll.il For 
this model, we find the distribution of the time before the fist collision tq and the distribution 
of the random normal Vi, which describes the first collision. Then we use these results to 
derive the simplified model, defined by ([I]) and (j2J). 

Since the medium is assumed to be isotropic, it is natural to suppose that the obstacles' 
centers form a Poisson point process, with the control measure proportional to Lebesgue 
measure A3. The only shortcoming of this assumption is that obstacles can intersect. We 
denote this Poisson process by II and write its control measure in the form (7t\r 2 )~ 1 \ 3 , where 
r is the sum of radii of the particle and an obstacle; A > is a parameter. We will see that 
A signifies the mean free path of the particle. 

In addition, we assume that for the initial speed vq the following condition holds: 

\v^\ 2 >r\a\, (3) 
where Vq is the projection of vq on the orthogonal complement of a. 

2.1. The distribution of To. We start with the following notations: for a v G R 3 , define 
the unit hemisphere 

S v := {u e R 3 : |w| = 1, (u,v) > 0}; 
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A tl ,t 2 ■= Y> s + — + rS V0+as , s G [t x ,h] 

Let us fix a t > 0. The inequality r > t is equivalent to the absence of obstacles' centers 
in the set Aoj, whence 

P{r >t}= P{n(A,t) = 0} = exp{-(7rAr 2 )- 1 A 3 (A, t )}- 
We claim that for volume of the set ^4 ,t ? which has the form of a curved "cylinder", 

A 3 (A,t) = vrr 2 f \v + as\ds (4) 
Jo 

(the naive explanation is the following: the factor of the integral is area of the "cylinder's" 
cross-section, and the integral is the "cylinder's" length). Then, for the distribution function 
of r , 

P{r >t}= e- x ^ & \^+ as \ ds . (5) 

Also, P{t < dt} = A _1 |fo|tit + o(dt), as dt — > 0, thus we see that the parameter A signifies 
the mean free path of the particle. 

Let us prove (j3J). Without loss of generality, assume that = (0, |^|, 0) T ; since 
a = (0,0, |a|) T , the first coordinate of the trajectory of the curve s i— > vos + as 2 /2 is zero. 
Denote by L{t) := L \vq + as\ds the length of this curve, and let s i— > (0, 72(5), 73(s)) T be the 
natural parametrization of the curve (that is a parametrization such that 72(s) 2 + 73(s) 2 = 1 
for any s). Then we can represent Aqj in the form 

Ao,t= {(0,72(s),73(s)) T + r«S( ,72( fl ),73( fl ))T, se [0,L(t))}, 
and, finally, introducing the function 

/ \ / 1 \ / sin^cos^ 

Q(s,<p,6):= \ 72 (s) +r i 2 (s) -73(a) cos^ 

\ 73(s) / \ 73 (s) 7 2 (s) / \ sin6»sinv9 

we have 

= g([0,L(t)),[0,27r),[0,7r/2]). 

By simple, but tedious calculations, it follows from (J3]) that for any / > the mapping 
Q : [0, /) x [0, 27r) x [0, 7r/2] — ► M 3 is bijective. Further, using the equalities 72(s) 2 + 7s(s) 2 = 1 
and 72(s)72(s) + 73(5)73(5) = (the second inequality is the derivative of the first one), the 
reader can easily prove that for the Jacobian of Q it is true that Jac Q(s, if, 9) = r 2 sin 9 cos 9. 
We finish the proof of PJ integrating the Jacobian over the set [0, L{t)) x [0, 2ir) x [0,7r/2]. 

There exists a very convenient representation of r . Let rj be an exponential r.v. with 
mean A; recall that F : M. 3 x M + — > R is defined as the solution of the equation 



I 



F(v,t) 

\v + aslds = t. 
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Since F(v, •) is monotone and F (v, J* * \v + as\ds^ = t, we have 

T = F(v ,v)- (6) 

2.2. The distribution of v\. Recall that at a collision the particle's speed V\ = v + ar 
changes to V\ — (1 + «)(Vi • vi)i>i, where v\ is the inner unit normal to the first obstacle 
at the point of collision. To simplify the notations, let us write v instead of v\. As in the 
previous subsection, we assume that v£ = (0, |«o-|,0) T . 

At first notice that v G Sy x \ this vector could be defined by its spherical coordinates 
(Vi/,^ 1 )) where the longitude (p v G [0, 2ir) is the angle between (1,0, 0) T and v, and the 
latitude 9^ G [0, ir/2] is the angle between V\ and v\. For any < t x < t 2 , ip G [0, 27r], and 
9 G [0,tt/2], define 



Ai.ta.v.fl := {v s + ^-+r{u G : ^ u < <p, 9 v u 0+as <9},sE [ti,t 2 )| C A tl 

Then 

p{</?„ < </?, ^f 1 <#,r G [t,t + dt)\ 

i^<^<0\r o e[t,t + d t)} = -L p{roe[M + , t)} i 

P{U{A ,t) = 0, UiA^t+dt^e) = 1 } 



*2- 



p{n(A, t ) = o,n(A,^) = i} 



+ o(l), rft^O 



(the term o(l) appears, because since A 3 (^4f t t + dt) = 0(dt), we have P{n(„4 tit+dt ) > 2} = 
o(dt)). The mapping Q from the previous subsection is bijective, therefore the sets Ao,t and 
At,t+dt are disjoint. Consequently, 

FU < <p,0? <0 r = t}= hm P (n(A^) = l} = X 3 (A^) 

but the numerator is the integral of | Jac Q | over |X(t), L(£ + rft)) x [0, </?) x [0, 9], the denom- 
inator is the integral of | JacQ| over [L(t), L(t + di)) x [0, 2-n) x [0, 7r/2], and we get 

P{^<^<e|v„ = i} = ^. 

We see that the distribution of v\ is invariant under rotations around Vi, and 

P{^i < 9} = sin 2 0; 

moreover, 6^ and r are independent. 

Let us find a suitable representation of v\. Suppose that a is uniformly distributed on the 
unit sphere S 2 and is independent of V\ (and of To); let v be the unit vector directed along the 
bisectrix of the angle between a and V\. Then, for any fixed Vi, the conditional distributions 
of v\ and u coincide. Indeed, both of them are invariant under rotations around Vi, and 
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P{0^ < 0} = P{0^ < 20} = (1 - cos 20) /2 = sin 2 0. Finally, from v)v = \{V t + \V x \<t) 
it follows that 

Vl - (i + a )( yi . ^ 1 y x _ i^(Vi + |Vi|<r). 

2.3. The simplified model. We simplify the model by stating that we can apply the results 
of two previous subsections to describe the particle's motion after each collision (i.e., we can 
simply replace v by the speed after the collision). Thus we obtain the model defined by 
and @. 

The main disparity between the models was discussed in Subsection 11.21 In addition, 
the models differ, because the distributions of tq and V\ were derived under the technical 
condition (j3J); clearly, the variables V n — ^^-(V n + \V n \a n ) and v do not have to satisfy it. 

However, our simplifications look quite natural. Indeed, the "disappearance" of obstacles 
after collisions in a sense means that a particle never returns to already met obstacles. This 
is reasonable if the obstacles are rare, because there is a drift in the direction of a. It is also 
sensible to neglect if r or \a\ is small; note that since distribution of V n converges to 7Ty, 

the variables (V n — ^-^{V n + \V n \a n ))~ L converge to a nondegenerate limit. 

3. Starting the proof of Theorem CD 

3.1. The position of the particle at the nth collision. In this subsection we shall find 
a suitable representation for X n := X(t n ), where t n := Yl^Zo T i ls the moment of the nth 
collision; additionally, put to := 0. 
Consider the new Markov chain 

*.-.= ( v "). »«; •»==( "° . 

\VnJ \-Va/\v \ 

It will be obvious later that the current initial condition describes the particle with initial 
speed Vq (at time zero happens the dummy collision, which does not change the speed). Note 
that V n and cr n are independent. Further, introducing the notations 

x:=v— ^—^-(v + \v\a), x = ( V ) e X : = M 3 x S 2 , 
2 \a/ 

we rewrite and (jSJ) as 

K+i = K + aF{$ n ,r ]n ). 
Thus for the new chain it is true that 

f& n + aF($r ni ri n )\ 
$ n+ i = , n > 0. 

\ &n+l J 

Let us agree to write coordinates of vectors of M. 3 using superscripts; recall that a 3 = \a\. 
From the trivial equalities V^ +1 = $ n , V^ +1 = $ n , and V 3 +1 = $„ + \a\r n it follows that 



X n +i — X n + 



( $ T \ 

\ ((Kf + i) 2 -(C 3 ) 2 )/(2|a|) / 



Xn + ]a~\ V n+l V n+l ~ Jn V 

\ ((K 3 + i) 2 - (CV)/2 / 



V. VYSOTSKY 



whence 



X 



n+l 



Besides, 



1 

\a\ 



n 

E 

i=0 



2 -- _3 , ^2- -3 

$ n+ l - $ $0 

V ((CI 3 ) 2 - ($; 3 ) 2 )/2 y 



+ 



1 n+l 

nE 



/ 



i=l 



^1^3 
^2^3 

V ((^ 3 ) 2 - (^ 3 ) 2 )/2 J 



(7) 



Ti 



1 v — v o -^3 1 /- —-3 ^~--3 x 1 



i=0 



n+l 

E 

i=l 



(8) 



Finally, denoting 



1 



1 3 ^3 
V V — XX 



/(*) : = ITT 



v 2 v 3 



^-2^3 
XX 



from ([7]) and ([8]) we get 

/ 

l 

X n — C\_at n 



(x 3 ) 2 )/2 














x 




i f 3 — x 3 y 





ex, 



l 

lal 



— -1- — 3 —1—3 
$ n $ n - $ $0 
. — 2- — 3 —2- — 3 
$„ - $ $0 



\ 



n 

V ((C 3 ) 2 - (%) 2 )/2 J 



Cl 

— -, — ral 



j=l 



($„ -$ )+^[/-c 1 / l ]($ i ). 
By definition, put g := / — c\h (the value of ci will be defined later) 



(9) 



3.2. The problem in terms of the Markov chain $ n . Between collisions the particle 
moves with constant acceleration, thus the process X(t) is defined by its values X n at the 
points t n . In fact, we can find the values of X l (t) and X 2 (t) by linear interpolation and the 
values of X 3 (t) by quadratic interpolation with the leading coefficient \a\/2. Analogously, the 
process Y t (s) is defined by its values at the points t n /t, but for E t 3 (s) we shall use quadratic 
interpolation with the leading coefficient |a|t 3//2 /2. 

Let Y t (s) be the following process: at the points t n /t put 



Y 



t 



X n — C\at r 



n > 



and define the values at other points via linear interpolation. The only difference between 
Y t (s) and Y t (s) is in the method of interpolation for the third coordinate. It is easy to 
see that Y t 3 (s) - Y t 3 (s) = |a|t 3/2 (s - t n )(s - t n+1 )/2 for s e [t n ,t n+1 \. Thus, denoting by 
n(t) the (random) number of collisions by the time t, for the norm || ■ ||c of the space 
C[0,1] =C([0,1],M 3 ) we have 



\Y t (-)-Y t {-)\\ c < max sup 

0<k<n(t) t k <s<t k+1 



\Y(s)-Y(s) 



max r k . 



%\ft 0<fc<n(t) 

Then, we introduce the process Z t (s), putting at the points t n /t 



* i=i 



(10) 
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and denning the values at other points via linear interpolation. Trajectories of Z t (s) and Y t (s) 
are piecewise linear and their points of interpolation have the same x-coordinates (namely, 
t n /t), therefore from (jUJ) we have 

\\Zt(-) - Y t (.)\\ c < -1 max |^(I<K | 2 + I^H + <*(\* k \ + |^|)| . (12) 

yt l<k<n(t)+l u J 

From (flOl) and (TT2T) we see that for proving Theorem [1] it is sufficient to check that for 
all initial conditions $o = % G X, it is true that 

— = max n, — >■ 0, — = max <P/J — > (J, r — >■ oo, (13) 

^/i Q<fc<n(t) ^/t l<k<n(t)+l 

and there exist constants c\ > and C2, C3 > such that in the space C[0, 1] 



z t (.)-^y(-), «-oo. (14) 

Thus we must study the properties of the Markov chain $ n in detail. The necessary 
facts from the Markov chain theory are stated in Section 01 In Section [5] we prove that the 
chain $ n possess some useful properties. These properties are used in Section [6] for proving 
([13]) and ([141) . The assertion (fHj) is the most difficult; we prove it applying the FCLT to the 
sequence g($> n )- 

The problem is reduced to studying properties of the Markov chain <3> n . 

4. Basic facts on Markov chains 

The purpose of this section is to describe conditions under which a Markov chain satisfies 
the law of large numbers (LLN) and the FCLT. We also give a simple method for checking 
this conditions. All the statements and definitions are taken from pjQ] ; in this section multiple 
references to this source are omitted. 

We begin with several notations. Consider a Markov chain $ n , with an arbitrary state 
space X equipped with a locally compact, separable, metrizable topology and Borel cr-field 
B(X). Let P(x, •) be the transition function of $ n , let P n (x, •) be the n-step transition 
function, and let ir be the invariant measure of the chain (in the considered situations, there 
exists a unique invariant probability measure). Calculating expectations and probabilities, 
we indicate the initial distribution of the chain, i.e., £(<&o), with subscripts. For example, 
P x .{$„ E A] imply that £($ ) = $c an d ^v®n imply that £(<J> ) = n. All the considered 
functions are assumed to be measurable. Finally, by P denote the transition operator; recall 
that by definition (Pf)(x) = f x f(y)P(x, dy), for any functional / : X — > R. 

4.1. Definitions. A Markov chain is called irreducible if there exists a nonzero measure \x 
on B(X) such that 

H(A)>0 F x {3n eN:f n GA}>0, x G X, A e B(X); 

any measure satisfying this condition is called irreducible measure of the chain. 

An irreducible chain is called aperiodic if there does not exist a d > 2 and there do not 
exist disjoint sets E\, . . . , E d G B(X) such that 

1) for all x G Ed, P(x, E\) = 1, and for all x G E iy P(x, E i+ i) = 1, i — 1, . . . , d — 1; 
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2) /i^X \^ jj i?^ = holds for every irreducible measure /i of the chain 

(we modified the definition from [TU] using Proposition 4.2.2 and Theorem 5.4.4). 

We say that a chain is (weak) Feller if the function P( • , A) is lower semicontinuous for 
any open set A C X. 

Let /i be a signed measure on £>(X), and let / : X — > [0, oo) be a functional. We define 
the f -norm of fi as 

:= sup / gdfi = / /d|/u| 
g-\g\<fJx Jx 

(the inequality |g| < / is pointwise). The 1-norm is called the total variation norm; the 
notation \\fi\\x is replaced by 

Let Pi and Pi be Markov transition functions, and let U : X — > [1, oo) be a functional. 
By definition, put 

hip pi,, \\Pi(x,-)-P2(x,-)\\u 

" i — -T2 [/ •— SUp — — -rr . 

xex U[x) 

A Markov chain <3> ra is ergodic if there exists a measure ir such that for any x G X 
it is true that ||P n (s,-) — 7r|| — > 0, as n — > oo; this yields that 7r is a unique invariant 
probability measure. A Markov chain is U -uniformly ergodic if there exists a measure 7r 
such that |||P n — n\\\u — ► 0, as n — > oo (we formally put vr(a;, •) := 7r(-)). Note that if a 
chain is [/-uniformly ergodic, then it is ell- uniformly ergodic, for any c > 1. For irreducible 
aperiodic chains the 1-uniform ergodicity is equivalent to the well-known Doeblin condition 
(see Theorem 16.2.3). 

Let g : X — > M. be such that g G -^ 1 (vr) = L}(X, B(X), n). The functional equation (in 
unknown a) 



g- Pg = 9 - / fi^vr (15) 

is called the Poisson equation. 
4.2. Theorems. 

Theorem 2. Lei $ n 6e an ergodic Markov chain, and let g G L 1 (7r). Taen /or any initial 
condition $ = a; G X 



1 /" 

lim Tj($i) = / odvr, P a 



-a.s. 



Proof. Follows from Theorem 17.0.1. □ 

Theorem 3. Let $ n be an irreducible, aperiodic, Feller Markov chain, and let Int(suppn) ^ 
for some irreducible measure /i o/$ n . Suppose that the Foster-Lyapunov condition holds: 
there exist a functional U : X — > [1, oo), a compact set C C X , and constants (3, b > sttca 

PC/(x) - < -^C/(x) + 61 c (a;), a; G X. (16) 
Taen $„ zs U -uniformly ergodic; moreover, U G L 1 (7r). 
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Proof. In view of Proposition 5.5.3 and Theorem 6.0.1, the first statement follows from 
Theorems 15.0.1 and 16.0.1; the last one is proved in Theorem 14.0.1. □ 

Theorem 4. Let $ n be aU -uniformly ergodic Markov chain, and let a functional g : X — > K. 
be such that g 2 < U . Then there exists a solution g of the Poisson equation (fl5l) : g G L 2 (n); 
and the constant 

il ■= [ (f - (Pg) 2 ) > (17) 

Jx 

is well defined. If j x gdn = 0, then for any initial condition $ — x £ X, in the space C[0, 1] 

St(s) := mm ± wMjw ^ ^ w{ . )t (18) 

where W is a Wiener process. 

Proof. The existence of a solution of the Poisson equation easily follows from Theorem 
17.4.2. The well-posedness of the definition of 7? (i.e., independence of the choice of the 
Poisson equation's solution g) follows from Proposition 17.4.1. The Cauchy-Bunyakovskii- 
Schwarz inequality implies that 7? > 0. Finally, the last statement is the combination of 
Theorems 17.4.4 and 17.5.4. Although in [10] the processes S t (s) are defined for positive 
integer t, in ([18]) the convergence over t G R simply follows from the convergence over 
t e N. ' □ 

5. Studying properties of the Markov chain $ n 
We shall frequently use the following trivial inequalities: for any x — ( v ) GXitis true 

that 

a\v\ < \x\ < \v\; 

recall that a G (0, 1). 

5.1. Irreducibility, aperiodicity, and the Feller property. First, let us prove that for 
any <3> = x the values of V2 = $1 + aF ($1,771) run over the whole 1R 3 . It is sufficient to 
prove that $1 runs through M 3 \ B(0, \x\) = [ti 6 I 3 : \u\ > \x\j, because for any fixed $1 

the values of F($>i,r]\) run through IR + . We certainly assume that rji and o~i run over the 
whole M. + and 5* 2 respectively. 

At a collision, the speed v G M 3 changes to v :— v — (1 + ce)(v , v)v (we temporary 
use the old representation). The reader will easily check that the inverse transformation is 
v = v — (1 + a~ 1 )({;, u)u, where v is the same as in the direct transformation. The values of 
v run through S v = {u G M 3 : \u\ = 1, (v,u) > 0}; consequently, we have (v,u) < 0, that 
is, v G <S_ij. Thus the speed after a collision could be equal to a v G M 3 iff the speed before 
this collision is contained in the set v^ 1 := {w G 1R 3 : u = v — (1 + a~ l )(v , u)u, v G S- v }. As 
above, define a as a unit vector such that v is directed along the bisectrix of the angle between 
cr and —v. Then a runs over S 2 , whence v^ 1 = [u G 1R 3 : u = v + — (\v\o~ — v), a G S 2 }. 

We shall prove that $1 runs over M 3 \ 5(0, |x|); recall that $1 = (^+ aF J^^)y Let us 
show that for any v G 1R 3 \ -8(0, \x\) the set f _1 + as, s > 0} is nonempty. This is 
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equivalent to the existence of a s > such that \x + as — — |— v\ = 1+ ° 1 \v\. But the 
left-hand side continuously depends on s and increases for large s. Since \v\ > the value 
of the left-hand side at s = is not greater than =i| — |i>|. Therefore the required s exists. 

Moreover, the distribution of $i has a density, and this density is positive on R 3 \ 
P(0, This fact is intuitively clear, because F(x,rj ) has a positive density on M + (see 
(J5J) and (jSJ)) and cri has a positive density on S 2 . The formal proof, whose main part is to 
calculate the Jacobian of the appropriate transformation, is omitted. 

Now it is obvious that the distribution of V 2 = $i + aF($i, rji) has a positive density on 
IR 3 . That is why the chain $ n is irreducible and \ 3 ®Usz is an irreducible measure; here A3 is 

the Lebesgue measure on M 3 and C/52 is the uniform distribution on S 2 . Indeed, $2 = (^2)' 

and the variables V2 and 02 are independent. 

We just proved that for any initial condition $ = x the distribution of $ 2) i- e -> -P 2 (a;, •)> 
has positive density with respect to A3 <8> Us*- Hence for any x G X the measures P 2 (x, •) 
and A3 ® C/52 are equivalent. Thus, by simple arguments, the chain $ n is aperiodic. 

To prove that is a Feller chain, it is sufficient to check that for any open set icK 3 
the function P( • , A X S* 2 ) is lower semicontinuous. Indeed, from (jSj) it follows that 

/>oo 

P(x,Ax5 2 )= / l A (x + at)X- 1 \x + at\e- x ~ 1 ^^ +a '\ d 'dt, 
Jo 

and applying the Fatou lemma, we obtain lower semicontinuity. 



5.2. [/-uniform ergodicity. 

Lemma 1. Let r\ be an exponential r.v. with mean A; then for any cGK 

sup Ee cF{v ' v) < 00. 

Corollary. For any initial condition $0 — % G X , the variables r n have exponential moments 
of any order. Moreover, 

supE x e CTn < 00. (19) 

n 

Proof of Corollary. Since r n = P($„,, r] n ) and $ n is independent of T] n , the proof is obvious. 

□ 

Proof of Lemma [I], We only consider the nontrivial case c > 0. Take an s > and a 
t> G M 3 . Since |— a\v \/\a\ + as\ = \\a\s — \v\ I < |i> + as|, we see that for any t > 

* . 

— ap|/|a| + as Ids < / \v + as\ds. 
io 
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Then, using the definition of F, we have F(v, •) < F(— a\v\/\a\, •) and therefore F{F(v,r}) > 
t} < IP^-F^( — a.|v|/|a|, 77) > t\. The right-hand side could be easily calculated, and 



poo 

Ee cF(v, v ) = _ e ct dF{F(v,7])>t} 
Jo 

POO 

= c e ct ¥{F(v,r]) > t}dt - 1 
Jo 

poo 

< c e ct ¥{F(-a\v\/\a\,7]) > t}dt - 1 
Jo 

POO 

= c / e d e- A " 1J, ° l|a| * -| " ll<fc tft-l 



|ti|/|a| poo 

e ct-A- 1 (|i;|t-|o|t 2 /2)^_|_ c / e ct-X- 1 (\a\t 2 /2-\v\t+\v\ 2 /\a-\) dt _i_ 

J\v\/\a\ 

We now estimate the integrals. For the first one, use the following inequality: if \v\ > 
4Ac, then ct - A -1 (|v|t - \a\t 2 /2) < -ct, for all t G [0, |u|/|o|]. Therefore if \v\ > 4Ac, then 

|i)|/|a| /"M/I a l 

e ct-X-HHt-\a\ty2) dt< / e ~ct dt< I e -ct dt<OQ 



and this bound does not depend on v. If \v\ < 4Ac, then 

r\v\/\a\ p4Xc/\a\ 

/ e ct-X-H\v\t-\a\tV2) dt < / e ct+X-l\a\ty2 dt < ^ 

Jo Jo 

Estimating the second integral, we apply the following: if \v\ > 5Ac, then ct— A _1 (|a|t 2 /2- 
\v\t+ |f| 2 /|a|) < — ct. Thus if \v\ > 5Ac, then 

OO POO POO 

e ct-X-H\a\ty2-\v\t+\v\y\a\) dt < / g-^f < / e ~ Ct dt < OO. 
\v\/\a\ J\v\/\a\ 

If |f I < 5 Ac, then we have 

e ct-X-H\a\t*/2-\v\t + \v\y\a\) dt< , ( , x . 



\v\/\a\ 



□ 



To prove the [/-uniform ergodicity of $ n , we apply Theorem [31 Take a c > and check 
that the Foster-Lyapunov condition ffl6|) holds for := e c ' x '. 
For the transition operator, 



PU{x) = E x U{®i) = E-^ 1 * 11 = E x e c|yi 2 (H+IHMI. 
We define the function 

T (H) := / e^-^^+MOl-^l^^^ t> g K 3 , 
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which is obviously monotone and j(\v |) — > as |t>| — ► oo. Because o"i is independent of V\ 
and of $ ) 

PU(x) = E, T (|T4|)e c l^ 

< E e 7(|Vi|)e c| *° l e c|a|7D 

< E x (7(l^/2|)l {|aN <^ |/2} + 7(0)l {|a|TO>|ro|/2} ) e^V'^ 

< ( 7 (|J/2|)Ee c l a l F ^) + Ee c l a l F ^°)l { | a | F( ^ o)> | J | /2} ) e c ^ 

< ( 7 (aM/2)Ee c l a l^°) + E e c W^°)l {WF( ^ o)>aH/2} ) C/(x), 

where as usual x = (^). It follows from Lemma [H that the factor of U(x) tends to zero as 
\v\ — > oo, whence for any (3 G (0, 1) there exists an R > such that 

PU{x) - U(x) < -PU(x), xiC R := 5(0, R) x S 2 . 

Clearly, for some b > 

PU(x) - U(x) < -PU(x) + blc R (x), xeX. 

Thus the condition ffl6l) holds and, consequently, for any c > the Markov chain $ n is 
e c M-uniformly ergodic. 

5.3. The invariant measure. By definition of LZ-uniform ergodicity, there exists a unique 
invariant measure ir of the chain. Since for every n the measure P n (x, •) is a product of some 
measure on i3(R 3 ) and Ug2, for the limit we also have 

7T = 7Ty <g U S 2, 

where iiy is a probability measure on £>(R 3 ). By the reasons of symmetry, try is invariant 
under rotations around the third coordinate axis. 

Further, we claim that the measure 7iy has a density. Indeed, in Subsection 15.11 we 
proved that P 2 (x, ■) has a density; moreover, it could be shown that for all n > 2 the 
measures P n (x, •) have densities, that is, P n (x, •) -< A 3 (g [752. Thus, passing to the limit, 

7T = TTy (g C/52 -< A3 (g C/52 and 7Ty -< A3. 

Theorem [3] implies that e c ' x \ £ L 1 {tx) for any c, and we can prove the following 

Proposition 1. For any initial condition $>o = x £ X , the variables V n have exponential 
moments of any order. Moreover, 

supE^e^" 1 < oo. (20) 

n 

Proof. We begin with \V n \ = |$„_i + ar n _i| < |K-i| + |a|r n _i < • • • < \x\ + \a\ Yh=o T i- 
Then, via the Holder inequality and Lemma [1], 

71— 1 71— 1 

E x e clv " 1 < e c ^E x T\ e claln < e c|3?l TT(E x .e c|a|n7 " ! ) 1/n < e c|J| sup Ee^ 1 ™^'^ < oo, 
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hence the exponential moments exist. To prove fl20|) . we combine the trivial inequality 
e c|«] < e ca 1 \y\^ £ or y _ fu\ £ x, and the definition of e ca 1 '^'-uniform ergodicity of $ n : 

lim E x e c|Kl1 = lim I e c|u| P n (x, dy) = / e cH c/7r(y) < / e co " 1| ^d7r(y) < oo. 

n-»oc n->ooJ x J x J x 



□ 



6. Finishing the proof of Theorem CD 
Recall that we must prove ([TBI and (fT4l) . Let us start with the following lemmas. 
Lemma 2. For any initial condition $o = x G X, 

r n = 0(logn), |K| = 0(logn), P x -a.s. 
Proof. Applying the Chebyshev inequality and then using ( fl9l) . we have 



E e Tn 1 
F x {r n > 21ogn} < < — supE x .e CT ". 

n n z n 

Thus the first statement immediately follows from the Borel-Cantelli lemma. Similarly, we 
prove the second statement via ( 1201) . □ 



Lemma 3. There exists a c 4 > such that for any initial condition Q = x <E X , 

lim — = c 4 , ¥ x -a.s. 

n— >oo n 

Proof. Recalling (jHJ) and the introduced notations, we see that it is sufficient to prove 

— 3 

lim — — = 0, F x -a.s.; 

n^oo n 

the existence of a c 4 such that 

1 n 

lim — — r /i 3 ($ l ) = c 4 , Pjr-a.s. 



i=l 



and positiveness of c 4 . 

Since |$ n | < from Lemma [2] we immediately obtain the first statement. Further, 
we can apply Theorem [2] to prove the second statement, because |/i 3 (x)| = \v 3 — x 3 | < (1 + 
< (1 + oT^e^l G L 1 ^) and thus h 3 G -k^ 71 ")- By definition, put c 4 := |a| _1 J x h 3 d^. 
The proof of positiveness of c 4 , which is quite simple, could be found in [13]. □ 

6.1. Proof of ( 1X31) . For the r.v. n(£), which denotes the number of collisions by the time t, 
it is true that 

lim n(t) = 00, F x -a.s. (21) 

t^oo 

To prove this, assume the converse. Then we can find a k > such that with nonzero 
probability the particle collides with obstacles only k times. Thus the probability of t& = 00 
is nonzero that contradicts with the existence of exponential moments. 
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By Lemma [3] and fl2Tj) , 

lim = c 4 , F x -a.s., 
t^oo n[t) 

but since t n (g <t< £ n (*)+i — tn(t) + r n(t) , from Lemma [2] we have 

lim ^tt = c 4 , K-a.s. (22) 

Applying Lemma [2] once again (recall that < |14|), we prove ({TBI (to be precise, we 
prove Pjr-a.s. convergence, which is much stronger). 

6.2. Definition of the constants ci, 02,03. We put 

Ci := c^ 1 ^! -1 / / 3 c?7r; 

for the quite tedious proof of positiveness of ci, see [H]. 

Now we can easily check that for g — f — c\h it is true that j x gdn = 0. In fact, 
for the first and the second coordinates, this follows from simple calculations, where the 
representation 7r = iry <g> C/52 and symmetry of ity are used. For the third coordinate, we 
apply the equality c 4 = | a, | 1 f x h 3 dn. 

Let us define ci and C3. By \g[x)\ < 2.5|a|~ 1 (l + a^ 1 )!^ 2 + ci(l + a -1 )|x|, there exists a 
c > 1 such that |g(a;)| 2 < ce'^L Thus (the chain <3> n is ce' x ' -uniformly ergodic) the functionals 
g 1 , g 2 , and g 3 satisfy the conditions of TheoremHl and there exist solutions g 1 , g 2 , g 3 G L 2 (tt) 
of the Poisson equations. Define 



<2 := \/c 4 1 7 9 2 2 , c 3 := W Vt^* 



,3^ 



and also 5- := (g 1 ^ 2 ^ 3 ) 

K:= [ (gf -(Pg)(Pg) T )d7r. 
Jx 



In [13] we proved (using the axial symmetry of ny) that the matrix K is diagonal and 

^ 91 = ^2> tllUS 

/ c 2 
K = cJ c 2 , 
V c 2 

6.3. Proof of ( TT4|) . In this subsection the following proposition plays the key role. 
Proposition 2. For the processes S t (s), defined in (TT81) . 

Proof. It is sufficient to show that for any nGl 3 

(S t (-),u) — > v /ci(y (•),«), t -> 00. 
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On the one hand, 

y^{Y(-),u) i sJc,({c 2 U^ + (C 2M 2)2 + ( C3U 3)2) ^(.) = ^(1^)^(.). 

On the other hand, the functional (gr, u) satisfies the conditions of Theorem HI thus 

(s t {s),u) = -= >y%,u) w v)> t^OO. 



And since, obviously, (g, u) = (g,u), from (TTT1) we get 

ll,u) = I {(g,u) 2 -(P(g,u)) 2 )d'K= [ ((g,u) 2 -(Pg,u) 2 )d'K = (Ku,u). 
Jx Jx 

□ 

Let us put S t (s) := S t (c~[ 1 s); then S^c^s) = yc^5 c -ij(s), and in the space C[0, 1] 
(and, moreover, in C[0, 1], for every / > 0) 

S t (-)^Y(-), t^oo. (23) 
Therefore we will prove (|T4l) if we show that 

\\Z t {-)-S t (-)\\ e ^0, t^oo. (24) 
For this purpose, introduce the process u t (s), putting at the points t n /t 

and defining the values at other points via linear interpolation. Using the definitions of 
Z t (s) and S t (s), i.e., (TTTj) and (|T8|) . we have Z t (t n /t) = S t (n/t). But S t (n/t) = SAc^n/t) = 
St lut (tn/t)) , whence at the points t n /t the equality Z t (t n /t) = S t {u t {t n /ty\ holds. However, 

Z t {s) = S t (u t (s)) 

is true for every s ! In fact, trajectories of Z t (s) and S t (u t (s)) are piecewise linear (for the 
last one, as a composition of piecewise linear functions) and their points of interpolation 
have the same ^-coordinates (namely, t n /t). As we saw before, the values at these points 
coincide. 

We see that Z t (s) is obtained from S t (s) by the random change of time. Suppose 
K(OI|c <2; then 



S t (-)\\c = \\S t (u t (.)) - S t (-)\\ c < u I (Ih(-) - id||c), 

^1(0,2] 



where u is the modulus of continuity, St|r 02 i is the restriction of S t (s) to [0,2]. Hence for 
any 1 > 5 > and e > 

F x {\\Z t (-)-S t (-)\\e>e} < F x {\\u t (-) - id|| c > 5} + ¥ x {u~ , (5) > e} 

bt \[0,2] 

< ¥ x {\\u t (-) - id\\ c > 5} + supP^a;- I (5)>e}. 



t>0 ^[0,2] 
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Let us proceed to the limit as t — > oo and then proceed to the limit as 5 — > 0. Now it is 
obvious that (1211) holds if for any 1 > 5 > 

hmP s {|K(-)-id|| c >5} = (25) 

and for any £ > 

lim sup F x {u~ I (5) > = 0. (26) 



<>0 I [0,2] 

At first we prove (125]) . Writing m*(s) in the explicit form, we have 



IK(0-id||c 
thus, by (1221), 



sup 

0<s<l 



c^n(st) st 



t 



lim || n £ ( 

t— >oo 



— • S 



T n(st) 

-id|| C 



< sup 

0<s<l 



C47l(st) — St 



0. 



F x -a.s., 



which is much stronger than 

It remains to check (126]) to complete the proof of Theorem [TJ The family of probability 
measures {F x o 5't|r 02 ]}t>o on B(C[0, 2]) is relatively weakly compact. This follows from (1231) 
and from P^-a.s. continuity of ||St||c[o,2] in t G [0, 00). The space C[0, 2] is a Polish space, 
thus the relatively weakly compact family {P,. o SA, nol }t>o is tight. By the well-known fact 
(see |2J) about tight families of probability measures on B(C[0,2\), for any e > 



lim supP x o SAr n9 Ap : u p (5) > e \ = lim supP x ^ u~\ (5) > e \ = 



i>0 



I [0,2] 
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